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Abstract. Among various definitions of quantum correlations, quantum discord has 
attracted considerable attention. To find analytical expression of quantum discord is an 
intractable task. Exact results are known only for very special states, namely, two-qubit 
X-shaped states. We present in this paper a geometric viewpoint, from which two-qubit 
quantum discord can be described clearly. The known results about X state discord 
are restated in the directly perceivable geometric language. As a consequence, the 
dynamics of classical correlations and quantum discord for an X state in the presence 
of decoherence is endowed with geometric interpretation. More importantly, we extend 
the geometric method to the case of more general states, for which numerical as well 
as analytical results about quantum discord have not been found yet. Based on the 
support of numerical computations, some conjectures are proposed to help us establish 
geometric picture. We find that the geometric picture for these states has intimate 
relationship with that for X states. Thereby in some cases analytical expressions of 
classical correlations and quantum discord can be obtained. 
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1. Introduction 

Correlation is the relationship between different things, and is a pervasive phenomena 
in nature. It is the way by which we learn the external world, and is the bridge on which 
we communicate each other and transmit information from this end to the other. 

In the classical world, correlations have been well studied from the viewpoint of 
information theory (see, for example, [1]). However "quantizing" classical information 
or correlation is definitely not an effortless work. The difficulties crop up in that quantum 
information, unlike the classical counterpart, is encoded in quantum states which may 
not be orthogonal and thus may not be distinguished unambiguously, and moreover, 
quantum systems can be correlated in ways inaccessible to classical objects. 

One of the prominent features of quantum correlation is entanglement. Entangled 
states cannot be prepared with the help of local operations and classical communication 
(LOCC) and thus they are nonclassical. Entanglement is indeed an important aspect 
of quantum correlation and is a prerequisite for many tasks of quantum information 
processing [2]. Nevertheless, entanglement is not the only aspect of quantum correlation, 
and the notion of quantum correlation is more general than entanglement. For example, 
there exists quantum nonlocality without entanglement [3, 4, 5]. 

Various approaches, other than through entanglement, have been proposed to 
study the correlations in composite quantum system. The first attempt at quantifying 
quantum contents of correlations is due to Zurek. The concept of quantum discord is 
proposed and develops to a measure of how non-classical the underlying correlation of 
two quantum systems is [6, 7, 8]. The important issue is the existence of quantum 
correlations beyond entanglement in separable states. Quantifying classical correlations 
in a bipartite quantum state and splitting the total correlation into a classical and 
a quantum part was presented by Henderson and Vedral [9]. Later the Henderson- 
Vedral (H-V) classical correlation was shown to have an operational meaning: The 
regularization of H-V classical correlation is just the maximal amount of common 
random bits obtained by one-way LOCC operations in excess of communication invested 
[10, 11]. In [12] Oppenheim et al presented an operational proposal, which comes from 
thermodynamical consideration, to quantify quantum correlations (see also [13]). By 
considering the amount of noise required to erase the correlation, Groisman et al gave 
an operational definition of the quantum, classical, and total amounts of correlations in 
a bipartite quantum state [14]. Recently Modi et al proposed a unified view of quantum 
and classical correlations [15]. 

The above considerations shed new light on the properties of the correlation 
incorporated in composite quantum system. Following the division of total correlation 
into classical and quantum part, many works have been devoted to study the roles 
played by different types of correlation in quantum processes, and reveal the relationship 
between them. These studies involve fuzzy measurement [16], mixed-state quantum 



Geometric picture of quantum discord for two-qubit quantum states 



3 



computation speedups [17, 18], broadcasting of quantum state [19], complete positivity 
of dynamics [20, 21], complementarity and monogamy relationship between classical and 
quantum correlations [22, 23, 24] and dynamics of discord [25, 26]. 

However, there is no effective method to calculate the exact results of quantum 
discord and other measures of quantumness analytically. Unlike the measure of 
entanglement, the new paradigms of quantumness of correlations are measurement 
oriented. What should be done in these paradigms is to extract information about 
system A by measuring another system B. Given a bipartite quantum system in the 
state p"^^ , when measuring system B gives the outcome k with probability pk, system A 
would be in some postmeasurement state p^. For a complete measurement on system 
B, the Pk and the pk are the members and probabilities of an ensemble of the local state 
of system A, that is, p^ = J2kPkPk- The accessible information about system A with 
respect to the particular measurement is given by 5'(p^) — J2kPkS{Pk)^ where S{p) is 
the von Neumann entropy of a state p. The major obstacle is to maximize the accessible 
information, or equivalently, to minimize the average entropy S = '^kPkS{pk)i over all 
possible complete measurements performed on system B. The explicit analytical results 
of quantum discord are known only for very special cases: Bell-diagonal states [27], 
X-shaped states [28] of two-qubit system, and Gaussian states of continuous variable 
systems [29]. 

Considering this problem, we propose in this paper a geometric method to describe 
the quantum discord of two-qubit quantum states. The geometric method is based on 
the idea of quantum steering ellipsoid which is defined in [30] . Quantum steering ellipsoid 
is such an ellipsoid in three-dimensional real space that each point in the interior 
or on the surface represents a postmeasurement state of one qubit when particular 
measurement has been performed on the other qubit. We denote the quantum steering 
ellipsoid by (£. The available postmeasurement states are constrained by the Or 
in other words, the decomposition of one local state, say p"^, can only be performed 
in the (£ (including the surface), namely, p^ = 'YlkPkPk for each p^ G We call a 
postmeasurement ensemble optimal if this ensemble can minimizes the average entropy. 
We also call the optimal ensemble as the optimal decomposition of the local state. It can 
be shown that the optimal ensemble can only be found on the surface of C This situation 
can be compared with the optimal signal ensembles studied in [31], where the output 
states of a noisy quantum channel are restricted in a convex set A, and the optimal 
signal ensemble {p^^*, Pfe''*} is such that the Holevo quality, x = S{p) — J2kPk^iPk) with 
each pk G A and p = '^kPkPk, reaches the maximum on {p^^ , p'k'^}- 

For two-qubit X states, there are only two candidates for the optimal ensemble. 
We call them equi-entropy decomposition and quasi-eigendecomposition respectively. 
The geometric picture of these two forms of decomposition is clear: equi-entropy 
decomposition corresponds to a horizontal line segment, while quasi-eigendecomposition 
to a vertical one. Then the known results for Bell-diagonal states and X states can be 
"seen" in this picture. Subsequently, we study the dynamics of classical correlations and 
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quantum discord in the presence decoherence. In [26], it has been shown that there is a 
sudden transition from classical to quantum decoherence regime for some Bell-diagonal 
states undergoing non-dissipative decoherence. We generalize this result to the case of 
general X states. The sudden transition can even be "seen" in the geometric picture. 

Not only can the geometric method be used to recover the known results, but also 
it should help us seek the possible analytical expressions of quantum discord for more 
general states. Following this line of thought, we consider a class of two-qubit states 
that have more complicated forms than X states. For these states, we can not give a 
thoroughly analytical procedure to derive the classical correlations or quantum discord. 
However, numerical computations give us the interesting results. We find that, just like 
the case of X states, there are only two possibilities as to the optimal postmeasurement 
ensemble. One is equi-entropy decomposition, and the other, although not the quasi- 
eigendecomposition, has intimate relation with the quasi-eigendecomposition. In the 
former case, we can write out analytical expressions for classical correlations and 
quantum discord, while for the latter further research is needed to characterize its 
property. We think that these phenomena revealed by numerical work should not be 
accidental coincidences. If these phenomena can be verified analytically, they will give 
us a geometric insight into the quantum discord. 

In Section 2 we give a brief overview of the concepts of classical correlations and 
quantum discord. In Section 3 we introduce a very useful tool, quantum steering 
ellipsoid, for evaluating quantum discord of two-qubit states. We present in Section 
4 the geometric picture to evaluate and describe the quantum discord of X states. For 
X states undergoing decoherence, the dynamics of classical correlations and quantum 
discord is studied and depicted geometrically in Section 5. More general states are 
considered in Section 6 and 7. Several conjectures and numerical tests are presented 
therein. Section 8 concludes. 

Throughout this paper, the logarithm has base 2. Numerical computations are 
performed by using Mathematica 8.0. 

2. Classical correlations and quantum discord 

Consider a bipartite quantum system composed of particle A and particle B, which are 
possessed by Alice and Bob respectively. The state of the whole system is described by 
a density matrix p^^. Total correlation between particle A and particle B is usually 
measured by the mutual information, that is. 



where p"^ and are local states of A and B respectively, p^(^) = TiB{A)ip'^^)- Mutual 
information quantifies the strength of the correlation. For product state p"^^ = p^<Sip^ , 
the entropy S^p"^^) is additive, namely, S{p^^) = 5'(p^) -|- S{p^), and it follows that 




(1) 
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the mutual information for any product state is zero. For maximally entangled state, 
such as ^ YltZo N)a ® "with d the dimension of the Hilbert space of subsystem 
A or B, the mutual information reaches its maximal value, 2 log2 d. It is shown that 
quantum mutual information is just the minimal rate of randomness that is required to 
completely erase all the correlations in p^^ [14]. 

The total correlation I can be spht into quantum part Q and classical part C, 
namely, X = Q+C. There are several ways to define the measure of classical correlations. 
Here we adopt the definition given by Henderson- Vedral [9], which quantifies the 
information gained about one subsystem from the measurement on the other. 

Suppose that Bob performs POVM measurements on his particle B. The set of 
POVM elements is denoted by = {Mk} with ^ and Efc ^fc = 1- The 
probability that outcome k is obtained is given by 

Pfc = Tr[p^^(l®Mfc)]. 

The postmeasurement state of particle A that corresponds to the outcome k is 

p^ = -Tr^[p^^(l®M,)]. (2) 

Pk 

Considering all POVM elements M^'s, the postmeasurement states of particle A are 
characterized by the ensemble {pk,Pk}- Note that Alice's local state p"^ remains 
unchanged, namely, p^ = YlkPkPk postmeasurement ensemble. Or in other 

words, Bob's POVM measurements induce a decomposition of Alice's local state p"^ 
into the ensemble {pk,Pk}- 

The information about particle A that is acquired by Bob's specific POVM Ai is 
given by 

s{p^)-J2p>^s{pt). 

k 

The dependence on the measurement procedure can be removed by maximization over 
all possible POVMs. The classical correlation is then defined as 

= max [S{p^) - J2P^S{pt)] = S{p^) - minJ^PkSipt), (3) 

M ^ — ^ M ^ — ^ 

k k 

where the maximization and minimization are taken over all of Bob's POVM 
measurements. The left-arrow over C indicates the situation that Bob performs 
measurement to acquire the information about Alice's system. Similarly, if Alice 
performs POVM, M = {Nj}, we can define the information gained about particle B 
by measuring particle A as 

= m^x [S{p^) - J^PjSipf)] = Sip"") - nun 5^p,S(pf ). 

j j 

Generally, ^ C^, meaning that the classical information is asymmetric. 
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It is natural to define quantum correlation as the difference between total correlation 
and classical correlation, namely, 

=I-C^ 

= min5^p,5(p^) + S{p^) - S{p^% (4) 

k 

Similarly for = X — C^. 

Quantum correlation Q is also called quantum discord. Quantum discord, which is 
originally defined as the the difference between two classically identical (but quantumly 
distinct) formulas that measure the amount of mutual information of a pair of quantum 
systems [7, 6], aims to capture all the quantum correlations, not limited to entanglement. 
There is a fundamental difference between entanglement and discord for mixed states, 
although they are equivalent for pure states. A typical example of this is the separable 
states with nonvanishing discord [8]. Other forms of definition of quantum discord can 
be found in [32]. 

To obtain quantum discord or classical correlation, one has to make considerable 
effort to minimize the average entropy S = ^kPkS{p^) over all possible measurements 
on particle B. In the next section, we introduce a useful tools, quantum steering ellipsoid, 
which will help us to establish a geometric picture about these concepts. 

3. Quantum steering ellipsoid 

First we express the states and POVM elements in Hilbert-Schmidt space. Let 
p^^ be a two-qubit state shared by Alice and Bob. It can be written as p^^ = 
i /3=o ® '^13 where (Tq is the 2x2 identity matrix, cTj {i = 1,2,3) are Pauli 
matrices, and Ra/3 = Tilp"^^ {aa<^cr p)] are all real numbers. We arrange the 16 coefficients 
Ra/3 into a 4 X 4 matrix R = {Rap)- Note that i?oo is just the trace of p^^ and equal to 
one. We write M^., one element of Bob's POVM, as = Yla=o^k,aCa- Similarly the 
state p^ (see (2)) can be expressed as = | J2lt=oyk,a<^a- 

Let's define two four-components vectors in the row form, = {yk,o, Vk) and 
Xfc = (xfc,o, Xk), where yk = {yk,i, yk,2, yk,^) and Xk = {xk,i, Xk,2, a^^.s)- Note that 
ykfi = 1 for all k and yk is the Bloch vector of p^. Direct calculation shows the following 
equation. 

3 

Pk Yfc = Xfc i?^, pfc = ^ i?o« Xk,a, (5) 

where the superscript T means matrix transpose. (5) provides the relationship between 
Bob's measurement and the corresponding components in Alice's ensemble. For 
entangled states, the matrix R is of full rank, and the vector is in a one-to-one 
correspondence to the vector x^. In the following, we will treat i? as a full rank matrix. 
Some states with singular R will be discussed in Section 6. 
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Also note that (5) imposes constraint on the vector y^: although can represent 
any projective measurement, the vector y^ can not be arbitrary. For example, y^ can 
not represent a pure state unless p^^ is a pure state. It is pointed out in [30] that the 
allowed y^ must satisfy 

yk[R-^vR-']iykf^O, (6) 

where = {R~^)'^ and rj = diag(l, — 1, — 1, — 1). In fact, (6) comes from the 

requirement that each is nonnegative. Noting that the 3-component vector yk = 
iVk,!, yk,2, yk;i) is the Bloch vector of p^, we can see that (6) describes a ellipsoidal 
region in three-dimensional real space. It means that, for each allowed to appear in 
Alice's ensemble {pkiPk^i corresponding Bloch vector is constrained within an 
ellipsoid (including the surface). The ellipsoid given by (6) is called "steering ellipsoid" 
in [30]. We denote it by (£. 

The steering ellipsoid renders concrete geometric picture when we do the 



minimization of the average entropy S . Each point belonging to C corresponds to some 
p^. By noting that the £ is convex and entropy function is concave, we see that the 
minimal value of S must be attained on the surface of C For any point on the surface 
of that is, for any vector y such that the equality in (6) holds, the corresponding 
vector x (see (5)), must satisfy x1 = %{ + x\ + x^. Such a vector x represents rank-one 
element of Bob's POVM that can be taken to be proportional to the one- dimensional 
projector 11, namely, M = mil. With the factor m absorbed into the probability p, 
we say that any point on the surface of C is induced by Bob's projective measurement. 
Therefore in order to obtain optimal ensemble of Alice's state. Bob need only to perform 
projective measurements. This fact has been pointed out by Hamieh et al in [33]. We 
give a geometric description here. 



4. Quantum discord of two-qubit X states 

For general two-qubit state, £ is too complicated to be dealt with. However for a 
specific class of states, called X states, we will show that the geometric picture is very 
clear. In various situations X states have been used to demonstrate significant quantum 
phenomena, for examples, entanglement sudden death or birth [34, 35], dynamics of 
quantum and classical correlations [25, 36], and sudden transition between classical and 
quantum decoherence [26]. The density matrix of a general two-qubit X state is 

fa 

_ h ve^' 
^~ ve-'" c 
y u e-'f" 

where a + b + c + d= l and u,v ^ 0. It is required that ^ ad and ^ be to assure 
the positivity of the density matrix. By local unitary operations, the off-diagonal entries 



u e 


d 



(7) 
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can be transformed to real ones. Since all correlations are invariant under local unitary 
operations, it suffices to consider X states with all the entries of density matrix being 
real. However we will remain at the form given by (7) for later references. 

The matrix R is given by 

1 



R 



\ 








a + b 
— c — d 



2u cos fi 
+2v cos u 

—2u sin /i 
—2v sin i> 









-2u sin fi 
+2v sin u 

-2u cos jj, 
+2v cos u 





a — b \ 

+c — d 




a — b 
-c + d J 



(8) 



When det{R) = I6{bc- ad){u'^ 



7^ 0, the inverse R ^ exists. From (6), we can write 



the equation of the ellipsoid that is, 



/2 



/2 

y2_ 

'■2 



p ~^ 77 ~^ 79. 



/2 

y3_ 



where y[ = yi cos (p — y2 sin 0, 1/2 = yi sin + 1/2 cos < 
and three major axes are given by 



y'^ = y?.- ^3 with 



h 
h 



u + V 



^{a + c){b + dy 
\ad — bc\ 



\u — v\ 



^{a + c){b + dY 
ab — cd 



(9) 

(/i + z/)/2, 
(10) 



{a + c){b + d)' {a + c){b + d)' 

In the coordinate frame {O' y[ y'2 y'^), the ellipsoid € takes the standard form, which 
comes from the form in (O yi y2 ys) frame by translating along y^ axis and rotating about 
?/3 axis. In the (O yi y2 ys) frame, the center of £ is at the point O' with coordinates 
(0, 0, Y3) (see figure 1). The Bloch vector of is given by 



(0, 0, a + b 



d), 



which is represented by a point A on the y^ axis. Point A can not be outside £ due to 
the fact that a + b — c — d & [Y^ — l^jY^ + l^]. 

We now proceed to find Alice's optimal ensemble {p/c,Pfc}°^*- When Bob performs 
complete projective measurements, Alice's ensemble has two components, i.e., k = 1,2, 
and both and are on the surface of 2. In the following, we demonstrate a geometric 
picture to describe how to obtain Alice's optimal ensemble and thereby the value of S^^^^. 

Imagine a class of planes that contain y'^ (or 7/3) axis and intersect the ellipsoid <B. 
Each section is an ellipse. By noting li ^ I2, the largest section must be in the 2/12/3 
plane. Then the largest ellipse (LE) takes the form 




Figure 1. Schematic plot of quantum steering ellipsoid £ of a two-qubit X state. In 
(O Ui t/2 2/3) frame, the center of (£, denoted by O', is on the j/3-axis. Point A represents 
the state of qubit A. 

The 1/3 axis intersects LE at two points G and H. Consider a line parallel to y[ axis and 
passing through point A. This line will intersect the LE at two points E and F (see figure 
2). Remember that point A stands for Alice's local state p^. Then can be expressed 
as either of the two forms of convex sum: p^ = pcPc + PhPh or p^ = PePe + PfPf, 
where pc is the state corresponding to the point G and the probability pc = AH/ GH, 
and similarly for other states and probabilities. These two forms of convex sum lead us 
to the average entropies, 

Sgh = PgS{pg) + PhS{ph), (12) 
Sef = PeS{pe) +PfS{pf) = S{pe) = S{pf). (13) 

We call (12) quasi-eigendecomposition, meaning that p^, pc and pn have the same 
eigenstate. We call (13) equi-entropy decomposition, meaning that p^ = pePe + PfPf 
with S{pe) = S{pf)- 

Now we state our main result. With Bob performing POVM measurement 
on his particle, the optimal postmeasurement ensemble of Alice's state is given by 
{PGjPg] Ph,Ph} or {pe = 1/2, pe; Pf = 1/2, p^}, and the minimal value of S is 

Sj^in = mm{SGH, Sef}- (14) 

It follows that the classical correlation and quantum discord are given respectively by 
C^(p) = S'(p^) - mm{SGH, See}, 
Qr{p) = m.m{SGH. Sef] + Sip"") - 5(p^^). 
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Figure 2. Schematic plot of the largest ellipse given by {y'l/h 



1. Two 



pairs of points, {E, F) and (G, H), are the only candidates for Alice's optimal ensemble 



which will give the minimal value of average entropy S . 



The result (14) can be derived by using the conclusion in [28]. As pointed out in 
[28], there are two candidates for Bob's measurements which will induce Alice's optimal 
ensemble. In our notations, these candidates are denoted by 4-component vectors x-t: (i) 
x± = (|; 0, 0, ±|); (ii) x-t = (|, ±| sin^, ±| cos 6*, 0), where the measurement parameter 
9 will be determined latter. 

For case (i), it follows from (5) that 
P+ = a + c, y+-- 



1, 0, 0, 



a 



c 



P- = b + d, y_ 



1, 0, 0, 



a + c 

b-d 

b + d 



It is easy to see that this case results in the two points G and H in figure 2. 

For case (ii), we have p+ = = 1/2 and y-t = (1, y±) with Bloch vectors y± given 

by 



/ ±2u sm{e - ^)±2v sm{e + v) \ 
=F2u cos(^^ — /i) ± 2f cos(^^ + v) 
a + b — c — d 



Because 
namely, S{p^ 



the entropy of the corresponding state is equal to each other, 
= S{p-). It follows that the average entropy is given by p+5'(p+) + 



P-S{pJ) = S{p^) = 5'(p_). To obtain classical correlation, we will maximize |^+| over 
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the parameter 6. In fact the maximal value of is attained when 6 = (vr + /i — z/)/2. 
In this situation, Bloch vectors y± is given by 

y± = 2{u + v) cos 0, =f2(m + f ) sin 0, a + h — c — , 

where = (/x + z^)/2. It is straightforward to check that y+ and y_ just correspond 
to the vector and in figure 2 respectively. Thus we have proved that the two 
pairs of points, (G, H) and {E, F), stand for the only two candidates for Ahce's optimal 
ensemble. Now the problem of finding quantum discord for two-qubit X states is reduced 
to a simple geometrical one. The only thing we have to take into account is the steering 
ellipsoid (£ and the largest ellipsoidal section. 



4.1. Bell-diagonal states 

To appreciate the geometric picture, let's consider a specific class of X states, i.e.. 
Bell-diagonal states. Although the quantum discord of Bell-diagonal states has been 
calculated explicitly in [27], we would like to provide a more concrete interpretation. 

For a Bell-diagonal state pbd given by 

Pbd = ^2^^mO"m®^m' 

where to = 1, 1 ± ts ^ \ti =i= and it is assumed that tit2^3 7^ 0, the steering ellipsoid 
<£ has the standard form, that is 

^1 y| _ . 

il l2 f-S 

Moreover Alice's local state, represented by point A, coincides with the origin point O. 
Without loss of generality, we assume ^ |t2| ^ l^sl- It is not difficult to see that the 
LE is given by 

= 1 

f2 ^ f2 
Ll {-2 

According to the previous analysis, if Bob performs two-element POVM measurement, 
the minimal value of 5* is attained at the pair of points {E, F) or the pair (G, H). Note 
that OE = OF = \ti\, OG = OH = \t2\ and ^ jtsl- It follows that S^^^^ = S{pe). 
Or more generally. 



Si^ = mm{h{\t^\),h{\t,\),h{\t,\)}, 



with the function h{x) defined by 



1+x 1+x 1-x 1-x 
n[x) = — log — — log — - — , (15) 



for X e [0, 1]. 
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5. Geometric picture of dynamics of quantum discord 

As an application of our result, let's consider the dynamics of quantum discord 
or classical correlation. Recently, this problem has received considerable attention 
[25, 26, 36]. It has been shown that for some Bell-diagonal states passing through 
phase damping channel, the classical correlation can be unaffected by decoherence. And 
more interestingly, the dynamics exhibits a sudden transition from classical to quantum 
decoherence regime [26]. We will show that in the geometric picture these phenomena 
can be "seen" clearly even for general X state. 

When each particle of a two-qubit quantum system undergoes the phase damping 
process, the evolution of the state is expressed as 

2 

Pad =J2(^^^ K.j)p{K, ® Kj)\ 

where Ki = diag(7, 1) and K2 = diag(A/l — 7^, 0) are Kraus operators representing 
phase damping channel and 7 = e~'"* with F the phase damping rate. Here we assume 
that qubit A and B endure the same noisy environment. At initial time t = the 
steering ellipsoid €(0) is given by (9). At time t > 0, the ellipsoid is transformed to 
which is expressed by 

/2 /2 /2 

Vi , 1/2 , ^ = 1 



ii'hy ii'h)' II 

That is, with 7 decreasing from 1 to 0, the radius of the ellipsoid along axis and 
that along 2/2 axis decrease continuously from li and I2 respectively to zero whereas the 
radius along y'^ axis (i.e., axis) remains the same. Since h ^ h, the LE is given by 
[1/1/(7^^1 )]^ + [yU^sf = 1 in the time evolution. The LE will shrink to y'^ axis, namely, 
points G and H remain fixed and points E and F approach gradually to y'^ axis. Also 
note that Alice's local state does not affected by phase damping and thus the point 
A is fixed. 

We now show that the dynamics of quantum discord and classical correlation 

C'^ can be demonstrated clearly in the geometric picture. To this end, it suffices to 



consider the minimal average entropy S^-^^. See figure 3. There are only two possibilities 



with respect to the initial value of 5*^^;^. One is that at t = Alice's optimal ensemble is 
determined by the points G and H, and then ^^^^^(t = 0) = Sqh- At time t, the points 
E and F move to Ef and Ft respectively, while points G and H remain unchanged. 
It follows from OEt < OF that SEtFt > Sef > Sgh- So in this case S^^^ is always 
given by Sgh and remains invariant. As a consequence, classical correlation does 
not change during the time evolution. 

The other case is that initially Alice's optimal ensemble is described by two points 
E and F, that is, S^^^{t = 0) = Sef < Sgh- Time evolution will make SEtFt larger 



Geometric picture of quantum discord for two-qubit quantum states 13 





G 










/At 


At 




f t\ \ 






0' 






H 












Figure 3. Visual interpretation of dynamics of classical correlation and quantum 
discord. 

continuously, until the evolution reaches the critical time, denoted by i, such that 
SEfFt = Sgh- For t > t, we have SEtPt > Sgh and then S^^^{t > t) = Sgh- This 
means that after the critical time i the classical correlation does not change any 
longer. In a word, the phenomena presented in [25, 26] also arise in general X states. 

It is noted that the discussion presented above is not limited to phase damping 
channel. In fact, the geometric picture applies to any quantum channel that preserve the 
X form of the state, such as all unital channels in the canonical form (e.g., Pauli channel), 
amplitude damping channel, etc. Classical correlations may not remain constant in these 
more general cases. 

6. States with singular R 

The geometric method presented in Section 3 and 4 is based on the 3-dimensional 
quantum steering ellipsoid, which requires a nonsingular coefficient matrix R. To extend 
this idea to the case of singular R, we consider in this section two classes of states. One 
is the class of X states with det(i?) = 0, which is in fact the supplement to the content 
of Section 4. The other is such a class of states coming from mixing two pure product 
states. For these states, we put forward a conjecture about the geometric description of 
the quantum discord. 

6.1. X states with det(-R) = 

It follows from the coefficient matrix R given by (8) that when ad — be = or u = v, the 
determinant of R vanishes. Recalling the ellipsoid given (9) and the parameters given 
by (10) and (11), we have the following cases. 

When ad — be = and u v, we see that h = and (the 1/3-component of f^) 



Geometric picture of quantum discord for two-qubit quantum states 



14 



is equal to Y3. Then the elhpsoid degenerates to the elhpse. In the {O' y'^y^y''^ frame, 
the equation of the elhpse is 

,2 ,2 

12 ^ 12 
n '■2 

And the position of Alice's local state happens to on the origin point O' . Since h > I2, 
the minimal value of S is given by 



I + I + 1-/1, 1-/1 
7^ log -1^ 7^ log 



When ad — be = and u = f 7^ 0, it follows that ^2 = ^3 = and = Y^. We have 
two points given by y[ = li and y[ = —li, which will determine the optimal ensemble of 
Alice's state. 

When ad — be ^ Q and m = f 7^ 0, we have ^2 = 0, and the ellipsoid degenerates to 
a ellipse in y[ 1/3-plane, that is, 

/2 /2 

^ = 1 
12 ^ 12 



In this case, S^^^ can be easily obtained with reference to figure 2. 

The last case is that m = f = 0. In this case, the density matrix p takes the diagonal 
form. If ac? — 6c 7^ 0, the state is classically correlated and the ellipsoid reduces to the 
two points G and H in figure 2, which determine the quantity S^^^. If ad — be = 0, the 
state is a trivial product state. 



6. 2. Mixture of two pure product states 

In this subsection, we will find the quantum and classical correlations in such states 
that can be written as 

p = X 1^1) {tPil ® IV'2) (V'2| + (1 - A) |0l) (0l| ® 102) (02| , 

where A G [0,1], and and {i = 1,2) are the states of particle A and B 
respectively. Since the correlations remain invariant under local unitary transformations, 
it suffices to consider the states with the following form. 

p = A |0) (0| ® |0) (0| + (1 - A) 1^) ® \v) (<^| , (16) 

where \ip) = cos a |0) + sin a |1), and \{p) = cos/3 |0) + sin/3 |1) with a, (3 G [0, |]. 

A special case of (16), where both \ip) and \ip) are set to be |+) = -^(|0) + |1)), is 
discussed in [9] (see also [33]), and numerical evaluation is performed to inquire about 
the classical correlation therein. Here we consider a more general case, and will put 
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forward a conjecture about the exact value of quantum discord or classical correlation 
from a geometric viewpoint. 

Suppose that Bob performs POVM measurement on his qubit to acquire 
information about Alice's qubit. As stated earlier, we need only consider projective 
measurements. Let Bob's measurement operators be M_|_ and M_, that is, 

1 ^ ^ 

M± = - ± xia.j, ± X2ay ± xacr^, (17) 

with x\ + x\ + x\ = 1 / 4:. Bob's measurement will give the result "+" with probability 
and the result "— " with probability p_. The corresponding postmeasurement states 
of qubit A are and p"^ respectively. 

Define two 4-component vector x+ and x_ as 



x± = ^) = ±a;i, ±3^2, ±a;3 



Denote by y± = (?/±,i, y±,2, l/i.s) the Bloch vectors of pj. respectively. Then from (5) 
we have, 

p± = ^ ±a;i(l - A)sin2/3±a;3[A + (1 - A)cos2/3], (18) 

y±,i = ~{\^^ ~ ^^^2" 

±a;i(l - A)sin2asin2/3 ±X3(1 - A) sin 2a cos 2/3 1, (19) 

y±,2 = 0, (20) 

l/±,3 = — |^[A + (l-A)cos2a] 

±xi(l - A) cos 2a sin 2/3 ±X3 [A + (1 - A) cos 2a cos 2/3] |. (21) 
We can see from these expressions that 

l/±,3 + Z/±,itana = 1. (22) 

(22) means that the two points corresponding to Bloch vector y+ and y_ are located 
on the line L that lies in yiUs plane and passes through the point (0, 1) with the slope 
— tan a. Note that the Bloch vector of is given by 



(^(1 - A) sin 2a, 0, A + (1 - A) cos 2a) . 



Then the point A, denoting the state p"^, is also on the line L. See figure 4. 

Before using this picture to find the quantum discord of the state given by (16), 
let's add some remarks. 
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Figure 4. Geometric picture about the states given by (16). Point E and F are 
on the hne L and denote the postmeasurement states of qubit A. Point A denotes 
the local state of qubit A. Line segment EF can slide between point (0, 1) and point 
(sin 2a, cos 2a). 

Roughly speaking, the set of postmeasurement states of qubit A (i.e., and p^) 
are restricted on the hne L. For convenience, we denote the two states by point E and 
F respectively in figure 4. Assume that E is on the left side of point A and F on the 
right side of A. Then it should be noted that E can not be located on the left side of 
point (0, 1), because the length of OE can not be larger than one. For the same reason, 
F can not be on the right side of point (sin 2a, cos 2a). Then the line segment EF, 
which represents the set of all available postmeasurement states of qubit A, slides along 
line L between point (0, 1) and point (sin 2a, cos 2a). This picture is somewhat different 
from that presented in Section 4, where the steering ellipsoid takes up a fixed region 
for given state and does not depend on the choice of the measurements performed by 
Bob on qubit B. But here the line segment EF is "moving", in the sense that both the 
length and the position of EF depend on Bob's measurements. 

Now we propose a conjecture about minimal average entropy S^-^^. Given a 
two-qubit state (16), Bob performs two-element POVM measurement on qubit B. If 
Alice's postmeasurement ensemble {pk, p^}k=+.- minimizes the average entropy S , 
then S{pi) = S{p^). 

In other words, OE = OF is the necessary condition which must be satisfied in 
order that the average entropy S takes the minimal value. To test this conjecture, we 
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select randomly 1.5 x 10^ two-qubit states with form given by (16), and for each state 
calculate numerically the value of S^^^ and the corresponding measurement parameters, 
namely, Xi with i = 1,2, 3. From (19), (20) and (21), we get the coordinates of point E 
and F and also the length of OE and OF. In figure 5, we plot the value of OF — OF 
for the 1.5 x 10^ states. We see that \0F — 0F\ <^ 0. Numerical results confirm our 
conjecture. 



4. )xlO"'' - 



2.0x10"' - 



; 0x10"" - 



■ 0x10"' - 



Figure 5. Numerical test of the conjecture that OE = OF is the necessary 
condition for {pe,Pf} to be the optimal ensemble of p^. 1.5 x 10^ random states 
are tested and the values of \OE — OF\ are calculated. The results show that almost 
\OE-OF\ < 10-s. 

If the conjecture is indeed true, we need only consider the situation that 5'(p^) = 
S{p^) or OF = OF. It follows that 

To obtain 5^^^^, we need only to maximize OF or OF under the condition that 
OF = OF. It is not a difficult work. Following this line of thought, we obtain the 
classical correlation and quantum discord of states (16) and plot the results in figures. 

In figure 6 and figure 7, we plot the classical correlation and quantum discord 
for the states given by (16) with A = 0.5 and A = 0.7 respectively. They are very 
similar to each other, but it should be noted that the plot of in figure 6 is symmetric 
with respect to the parameter /3, while it is not the case in 7. To see this, refer to figure 
8. 
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Classical correlation for X=\/2 Quantum discord for X=l/2 



Figure 6. Classical correlation and quantum discord of the states given by 
(16) with A = i. 

It can be seen that among all states given by (16) with fixed A, the one with maximal 
classical correlation is of the form 

A|0) (0|® |0) (0| + (1 - A) |1) (1| ® |1) (1|, 

and the one with maximal quantum discord is given by 

A |0) (0| ® |0) (0| + (1 - A) |1) (1| ® 1+) (+1 . 

It is not difficult to see that the choice of A = 1/2 renders the largest classical correlation 
and quantum discord. 

7. More general states 

It is desirable to apply the geometric picture for a wider class of states. This section is 
devoted to extending the discussion in Section 4 about X states to a more general case. 
Recall that in the geometric picture of X states Alice local state is located on the 
(or 2/3) axis, which is one of the symmetric axes of the ellipsoid (£ (see figure 1). How 
about the case that point A deviates from axis? We will discuss in this section this 
type of states. To begin with, let's see an example. 

The example we will consider comes from [37]. Suppose a two-qubit pure state is 
given by [38] 
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Figure 8. Quantum discord of states given by (16) witli a = tt/6 and A = 0.7. 
It is not symmetric with respect to /3 G [0, f ]. 
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where iV'o)^ = + l^)) \i^i)A ~ I I*-*) + 1 1-'-)- qubit A pass through a 

quantum channel, the Kraus operators of which are given by 

^i = |0)(0| + -^|l)(l|, A, = ^\0){1\. 

Then output state of the channel is 

2 

1=1 

Note that p^^ here is not an X state. With Bob performing measurement, the steering 
ellipsoid is given by 

yl + yl + 2{y3-iy = 1- (23) 

The Bloch vector of p"^ is 

The ellipsoid (23) is symmetric under rotation about axis. The form is similar to 
that for X states. However, the vector does not lie on axis (see figure 9). We can 
not obtain the S'^j^ analytically, but numerical evaluation reveals an interesting result: 
Considering any line passing through point A and intersecting the ellipsoid at two point 
E and F, the minimal value of S is reached at such E and F that OE = OF, or, 
^ 

where Se or Sp have analytical expressions. We can see that it is the equi-entropy 
decomposition. Geometric description is clearly demonstrated in figure 9. 



The numerical value of S^^^^ is equal to 0.2804, and classical correlation and 
quantum discord is given by 

= 0.0118, = 0.0338. 

The above example motivates us to consider more general two-qubit states. That 
is, the quantum steering ellipsoid can be translated along y^ axis and no rotation is 
allowed. The position of local state p"^ is restricted in yi y^ plane. 

We take into consider the states p"^^ with the following R matrix. 

/ 1 Si S3 \ 

n til ti3 

t22 
\ rs t3i ^33 J 
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A /..--''\ 




E \ 














Figure 9. The ellipsoid 2; given by (23) is cut by a plane that is parallel to yi j/2 plane 
and contains point A. The intersection is a circle. Passing through point A, any chord 
of the circle leads to the optimal ensemble of p"^. 



It is assumed that R is non-singular, namely, det(-R) 7^ 0. Obviously, the Bloch vector 
of is given by = (ri, 0, r^). Assuming that Si,ti3 7^ 0, we choice the parameters 
til and as 

, _ n - 53^13 , _ rir^Si - ritgi + S3ti3t3i 
til — 5 ^33 — 



Si 



S^t 



1''13 



Under these conditions, we construct the states p^^ randomly. The quantum steering 
ellipsoid C given by 

yf , yl , iys^zl^ = 1 



where 



rf{l - si) - 2riS3ti3 + {sj + sl)t 



13 



slil-sj-sl) 



■-22 



-.2 ' 



[rj{l - sj) - 2riS3ti3 + t^sj + sl)]{rsSi - 



slt'l^{l Sj 



.2\2 



^3^1^13 - ^1^3(^3-51 - ^31) - ^13^3l(-Sl + ^3 
Siti3(l -sj- sl) 



Consider a class of lines passing through point A and intersecting the at point 
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M and N. It follows that = PmPm + PnPn- Corresponding this decomposition of 



p , the average entropy is given by Sj^^^ = PmSm + PnSn- We will find the minimal 



S among all these lines. Numerical results can be classified into the following two 
categories. 

Class I — The optimal line is parallel to yi axis or y2 axis. Denote by E and F 
the intersection points with (B. The minimal S is given by = Se = Sp, meaning 
that it is an cqui-cntropy decomposition. Moreover, let A' be the point that is on the 
line segment EF and symmetric to point A about axis. Then any point A" between 
the point A and A' has the same minimal value of average entropy, namely, 

sZ^^Se^Sf for A"eAA'. 

Class E — The optimal line is not parallel to yi y2 plane. In this case, let's consider 
the point A with the coordinate (0, 0, r^). Point A in fact corresponds to the projection 
of the vector r onto yz axis. We find that the minimal value of the average entropy 5" 
is given by the two points G and if, which are the upper and lower apex of the ellipsoid 
C respectively, that is, 
^ 

'S'min = PgSg + PhSh- 

It is a quasi-eigendecomposition. 

j\ 

In the case of Class I, S^:^^ have analytical expression. For Class I, we only see that 
it has relationship with quasi-eigendecomposition. More effort is necessary to acquire 
further insight. 

8. Conclusion 

We present a geometric method as to how to describe and evaluate the minimal average 
entropy, which is the major obstacle in the computation of classical correlations and 
quantum discord. For two-qubit states, the available ensemble of postmeasurement 
states of qubit A, which comes from the measurements performed on qubit B, is 
restricted in the quantum steering ellipsoid. The optimal ensemble can only be found 
on the surface of the ellipsoid. 

For two-qubit X states, the geometric method provides a clear picture as well as 
exact results. We show that for X states the optimal decomposition is alternative: 
equi-entropy decomposition or quasi-eigendecomposition. In the geometric picture, 
equi-entropy decomposition corresponds to a horizontal line segment, while the quasi- 
eigendecomposition to a vertical one. When an X state passing through some quantum 
channels, the dynamics of classical correlation and quantum discord can be easily 
analyzed in the geometric picture. 

We extend the discussion about X states to a more general case by relaxing the 
requirement that the reduced density matrices are of diagonal form. Little is known 
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about the classical correlations or quantum discord of these states. We perform 
numerical computations. A consequence of the numerical results is the following 
interesting alternative: the optimal decomposition is either equi-entropy decomposition, 
or is intimately related to the quasi-eigendecomposition. Combining these numerical 
results with the exact results for X states, we think that the geometric method may 
be generalized analytically rather than numerically. Further thought in this direction 
might be worthwhile. 

The geometric viewpoint presented in this paper offers an alternative way to 
interpret and compute classical correlations and quantum discord. It is also useful 
in elucidating issues related to decoherence. The remaining problem is to verify the 
conclusions we draw by means of numerical method. If these conclusions are indeed 
true, they will be helpful to work out the exact results of quantum discord analytically. 

This work was supported by National Nature Science Foundation of China, the 
CAS, and the National Fundamental Research Program 2007CB925200. 
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